QCD REVIEW 

GIUSEPPE MARCHESINI 
Universitd di Milano-Bicocca 
and 

INFN, Sezione di Milano-Bicocca, Italy 

This review is focused on QCD theoretical issues and their phenomenological relevance specially for 
LHC. It is incomplete and mostly neglects the phenomenology of long distance physics 



1 Introduction 

In 1973, with the discovery of asymptotic 
freedom ^, QCD was at the frontier of par- 
ticle physics studies. Now, in 2006, QCD is 
strongly installed (returning) at the center of 
particle physics researches. This, not only be- 
cause of the abundance of data on jet emis- 
sion at HERA and Tevatron, but especially 
for the necessity of preparing tools for discov- 
eries at LHC. Indeed, events with large Et, 
as in the decay of new massive particles, are 
accompanied by intense hadron emission. So 
the interpretation of a new elementary pro- 
cess requires an accurate description of short 
distance hadron physics, and this is the do- 
main of perturbrative QCD. Therefore calcu- 
lations have been performed in recent years in 
order to produce accurate hard physics pre- 
dictions for LHC. 

The results of high order QCD studies go 
beyond their phenomenological importance, 
they are exposing various new features of 
Feynman graphs which may point to new 
general properties of quantum gauge theories. 
Moreover the possibility that QCD can be 
viewed as a solvable string theory is getting 
strength after the discovery of the AdS/CFT 
correspondence. In the last two or three 
years this fact encouraged studies of "solv- 
able" string theories leading to phenomcnol- 
gical results in QCD. 

QCD studies are developing in many di- 
rections and I can discuss only a limited num- 
ber of points. Before describing NLO re- 
sults with their relevance for phenomenol- 



ogy and for understanding new features of 
QCD in general, I describe new avenues of the 
phenomenological attempts to connect QCD 
with string theory and its impressive list of 
results. 

2 Phenomenology of QCD as a 
solvable string theory? 

String theory originated from the Veneziano 
model ^ for strong interactions even before 
the QCD era. The possible relation of QCD 
with string theory is based on various key ob- 
servations. First, QCD Feynman graphs can 
be embedded ^ on the topological expansion 
of string theory (sphere, torus, etc.). Second, 
a consistent string theory must have more 
than four dimensions. Recently, Polyakov 
and Maldacena ^ suggested a holographic 
correspondence between a four dimensional 
gauge theory and a string theory in higher di- 
mensions. That is, the gauge theory observ- 
ables correspond to the observables on a four 
dimensional boundary of the string theory in 
higher dimensions. This correspondence was 
explicitly shown by Maldacena for the large 
Nc four dimensional supersymmetric Yang- 
Mill theory with four supercharges and the 
ten dimensional string theory with AdSs x 
metric (AdS/CFT correspondence). How to 
go from SYM theory to QCD? This question 
is one of the major issues in recent string the- 
ory studies ®, one tries to move away from 
conformal symmetry and reduce the symme- 
tries to approach QCD. 

In the last two or three years this prob- 
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1cm has been approached also directly form 
the QCD side (AdS/QCD correspondence). 
The question is: what are the characteristics 
of a string theory that holographically would 
reproduce in the four dimensional boundary 
the key proprieties of QCD. Results in this di- 
rection arc quite abundant. There are many 
approaches with a common starting point: a 
string theory with a AdS modified metric. 
In five dimensions (the minimum allowed for 
a possibly consistent string theory) typically 
one considers 

with Xn the usual four coordinates and z a 

fifth coordinate. The AdSs modified warp- 
ing factor h{z) accounts for confinement. It 
breaks scale invariance introducing a scale re- 
lated to the string tension. The QCD ob- 
servables are found at the boundary z = 0. 
At this point one tries to solve the string 
theory (by semiclassical approximation) and 
then deduces the QCD obscrvables typically 
at large distance (the region consistent with 
the semiclassical approximation used to solve 
the string theory). The results obtained are 
so numerous that I can only list them: 

i) masses, decay width, effective cou- 
plings, chiral symmetry breaking parame- 
ters for scalar, pscudoscalar, vector and axial 
mesons There are few parameters to fit and 
the agreement with the experimental data is 
reasonably good; 

ii) confinement with linear Rcgge trajectories 
and radial excitation number * . In a simple 
model one finds m^^ = R~^{J + n); 

iii) heavy quark potential at large and short 
distances; temperature phase transition and 
deconfinement in agreement with lattice cal- 
culations ^; 

iv) energy loss and jet quenching in heavy ion 
collisions ^°; 

v) high energy scattering at fix angles and 
exclusive form factors 

vi) Pomeron and the BFKL anomalous di- 
mension 



This long hst of AdS/QCD studies 
speaks for itself. String theory not only 
addresses the question of unifying particle 
physics with gravity, but starts to produce 
quite a number of phenomenological QCD re- 
sults. Whether this road is a correct one and 
QCD can actually be "solved" by string the- 
ory will be found out in the future. 

3 QCD expansion as string theory 
in twistor space? 

In 1986 Park and Taylor discovered that 

the many page long gg gggg amplitude at 
tree level and in the maximal helicity viola- 
tion (MHV) configuration can be reduced to 
a small single line formula. The formula was 
generalized 1^,14 ^.^^^ j,.^g^ with any num- 
ber of gluons. The multi-gluon amplitude can 
be written (a, X,p = colour, helicity, momen- 
tum) as 

perm 

with Tr(- • •) the colour order factors and 

M{Xi\2 • ■ ■ Xn) the helicity and momentum 
ordered amplitudes. For the MHV configura- 
tion (all positive helicities = +1 but two 
Xi = Xj = —l) one has the single line formula 



^MHV 



(PiPj) 



(1---) = ^^P^ T (1) 

(P1P2) ■ ■ ■ {PnPl) 

with {piPj) = y/2piPje^'^^^ and a phase 
proportional to the relative transverse mo- 
mentum. This simple formula, although valid 
only for the MHV tree amplitudes, indi- 
cates that Feynman diagrams have structures 
which are surprisingly simple. One is then in- 
duced to find the origin of this simplicity. An 
attempt in this direction was made by Wit- 
ten who suggested that perturbative gauge 
theories could be viewed as a string theory in 
the twistor space. The ground for this pro- 
posal is simple. Using the 2-spinors (massless 
Dirac equation solutions) one has 

{pp')=u-{p)u+{p') , Pf_c = ^u±{p)-yfiU±{p). 
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This shows that, while the momentum con- 
servation S'^{^pi) depends on both u- and 
u+ (together with the corresponding conju- 
gate spinors u+ the amplitude M^^^ 
depends only on (together with its con- 
jugate one U-), see (1). It is then natu- 
ral to perform a "Fourier transform" in U- 
and in the "Fourier" space, the twistor space, 
one finds that the amplitude with momentum 
conservation is described by a line in which 
the gluons arc attached in an ordered way. 
The problem is then, how to describe tree 
amplitudes which are not in MHV configura- 
tions and how to go beyond tree level. The 
results in these two directions have been re- 
ported by S.Moch and Z.Berm ^"^ . 

The prescription on how to go beyond 
MHV at tree level have been simplified and 
generalized i^.ig consists of sewing to- 
gether MHV amplitudes to construct the 
non-MHV ones. To do this one needs to 
use off-shell MHV amplitudes which arc ob- 
tained by analytical continuation in the com- 
plex momentum space. Tree amplitudes have 
been constructed for various processes: mass- 
less fermions Higgs boson EW vector 
bosons Results at one loop are also ob- 
tained The analytical expression for the 
one-loop gg gggg amplitude in all hclicity 
configurations will soon be obtained while its 
numerical evaluation is available 



4 High order QCD results 

Tree level amplitudes involving many QCD 
partons and next-to-leading order (NLO) am- 
plitudes are very important for LHC studies. 
Crucial studies are the search for the Higgs 
meson and for all signals which could indi- 
cate a way to complete the Standard model. 
The relevant events involve large mass par- 
ticles accompanied by intense hadron emis- 
sion. Thus it is crucial, for the interpretation 
of the events, to have accurate QCD predic- 
tions. This requires calculations of a large 
variety of many parton NLO matrix elements 



and distributions needed for: i) direct studies 
of new physics signals; ii) merging exact ma- 
trix elements with QCD resummation results 
for jet shape distributions and with Monte 
Carlo simulations of QCD jet emission . 

There are various numerical programs 
to compute many-leg amplitudes at tree 
level The most common tt decay {tt 
hhW^W~ — > bbqqq) involves 6 final state 
jets. Various processes have been computed 
at NLO and some presented at this confer- 
ence. The LHC "priority" wishfist for NLO 
calculations includes the following processes 



V e {z, w, 7} 


background to 


pp F jet 
pp — » ttbh 

pp ^ tt + 2 jots 
pp^VVbb 
pp^ VV-f 2jcts 
pp -\- 3jets 
pp^VVV 


ttH 
ttH 

tin 

VBF^ H VV, ttH 
VBF^ H ^VV 
new physics signatures 
SUSY trilepton 



The procedure for NLO calculation of an in- 
clusive distribution is in principle "simple". 
Typically, one starts from the tree level am- 
plitude for the process under consideration 
with n-partons (Born approximation). Then 
one computes the (n-l-l)-tree amplitude (real 
contribution) and the one-loop correction to 
the Born amplitude (virtual correction). Fi- 
nally one puts together both contributions 
to construct the distribution and check that 
collinear and infrared divergences cancel (for 
regular observables) . Since in general the 
inclusive sum is done numerically, the can- 
cellation needs to be controlled analytically 
first, a difficult issue which requires the un- 
derstanding of the physics of the problem 
(see for instance ■^^). Results and discus- 
sion on this issues have been presented at 
this conference Various techniques are 
used. On one hand there are semimimerical 
approaches. An example is the gg gggg 
one loop amplitude (some of the helicity 
configurations have been computed analyti- 
cally by twistor techniques). On the other 
hand there are direct analytical approaches. 
Powerful methods ^"^ are based on Cutkosky 
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rules (unitarity), on the use of Passarino- 
Veltman reduction of any one-loop ampli- 
tudes in terms of a basis of scalar integrals 
and on recurrence relations. 

The results of these studies are usually 
obtained after profound understanding of the 
general structure of Feynman diagrams. Of- 
ten one find general properties in gauge the- 
ories which points toward simple structures. 

Contributions to this conferences on phe- 
nomenological studies at Lep, Hera, Tevatron 
and LHC of high order results results have 
been presented. They are: 

i) running coupling measurements at Hera 
and Lep. Lattice "^^ calculations can be used 
to reduce the theoretical errors; 

ii) jet emission studies at Lep, Hera, Teva- 
tron and LHC. In particular at the Tevatron 
kt jet-finding algorithms start to be used; 

ii) heavy flavour production at Hera and 
Tevatron; 

iii) W/Z and W/Z+jcts at Tevatron ^5; 

iv) Higgs and W/Z production to NNLO ^6 
at LHC. 

The general comment for these analyses 
is that NLO corrections improve the accuracy 
and the description of the data. 

5 High order parton splitting 

Parton density functions enter DIS and, due 
to QCD factorization, hadron collider distri- 
butions. Fragmentation functions, which de- 
scribe inclusive final state emission, enter all 
collider studies. Their Q^-evolution is gov- 
erned by the corresponding space- and time- 
like parton splitting functions (anomalous di- 
mensions) which have been computed in 
1980 at two loops both for the singlet and 
non-singlet cases. Recently the anomalous di- 
mensions have been computed at three loops: 
the singlet and non-singlet ones for the space- 
like case ^^; the non-singlet ones for the time- 
like case These very important results ob- 
tained in MS scheme have been already used 
for various phenomenological studies ^'': Su- 



dakov resummations, Icpton pair and Higgs 
boson production, quark form factor, thresh- 
old resummation, DIS by photon exchange, 
longitudinal structure function, non-singlet 
analysis of deep inelastic world data. 

High order anomalous dimensions are 
also important to understand general features 
of QCD and gauge theories in general. I 
discuss here two examples: relation between 
DGLAP and BFLK evolution in SYM the- 
ory and relation between space- and time-like 
anomalous dimensions. 

DGLAP and BFLK evolutions. It 
has been shown that in the AT = 4 Super- 
symmetric Yang-Mills theory there is a deep 
relation between the BFKL and DCLAP evo- 
lution equations. In this theory, the eigen- 
values of the space-like anomalous dimension 
matrix are expressed in terms of a univer- 
sal function constrained (obtained) from the 
BFKL equation. This was checked at two 
loops by direct calculations and at three loop 
with the anomalous dimensions obtained for 
SYM from the QCD ones It is important 
to explore to what extent the relations be- 
tween BLKF and DGLAP can be extended 
to QCD. 

Relating S- and T-evolution. The 

search for a relation between S- and T- 
anomalous dimension (space-like j-{N) and 
time-like 7+ (TV)) has a long story: Drell-Levi- 
Yan relation Gribov-Lipatov relation 
the analytical continuation 37,44_ Consider 
DIS with q the large space-like momentum 
transferred from the incident lepton to the 
target nuclcon P and e+e~ annihilation with 
q the time-like total incoming momentum and 
P the final observed hadron. The Bjorken 
and Feynman variables in DIS and e+e~ are 
-2 2{Pq) 



Xb 



Xf 



2{Pqy r 
These variables are mutually reciprocal: af- 
ter the crossing operation P — > —P one x 
becomes the inverse of the other (although 
in both channels < a; < 1 thus requiring 
the analytical continuation). This fact was 
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the basis for the seareh of reciproeity rela- 
tions between 7-(VV) and 7+(VV) (here N is 
the MeUin moment conjugate to x). 

Recently it has been noticed that new in- 
formation on the relation between 7^ (TV) and 
7+(-/V) could be obtained by taking into ac- 
count that such a reciprocity property x 
1/x (x — xb or xp) can be extended to the 
Feynman diagram for the two processes and, 
in particular, to the contributions from mass- 
singularities described by multi-parton split- 
ting. Consider the three-parton vertex kine- 
matics of the decay fco in the DIS sit- 
uation: ko <0,k'^ < 0, fc'^ > 0. To change to 
the annihilation kinematics, —k — > —fco -I- fc', 
one has to change signs of fc+o and fc-j. and 
of the corresponding virtualities, fcg > 0, 

> 0. The virtuality fc^ enters the de- 
nominators of the Feynman diagrams. In or- 
der for the transverse momentum integration 
produce a logarithmic enhancement, the con- 
ditions must be satisfied 



\kl\<^\k'\ = z'^.n\ 



a = ± (2) 



with a = —l and cr = 1 for S- and T-case re- 
spectively. This kinematical fact has a strong 
impact on the relation between the S-and T- 
probability Da{N,K^) to find a parton with 
virtuality up to k^. Prom (2) one directly 
deduces the following reciprocity respecting 
equation (RRE) 

d^2DAN,^) =j^{N)D^{N,K^) (3) 



= / —z^P{z,as)D4N,z''K 
Jo ^ 

The difference between the two channels is 
simply in the fact that the virtuality of the in- 
tegrated parton distribution is z'^ , see (2). 
The splitting function P{z, as) does not de- 
pend on the S- or T-channel (its Mellin mo- 
ments arc not the anomalous dimensions). 
The running coupling in the splitting func- 
tion depends on the virtuality in a reci- 
procity respecting form. This equation (in 
general a matrix equation) is non-local: for 
a — — I (ct = 1) the right hand side involves 



^ dz 



the parton distribution with all virtualities 
larger (smaller) than k'^ . So RRE is not suit- 
able for explicit calculations of the anomalous 
dimensions, but for relating them. 

Is eq. (3) correct? It is the result of the 
vertex kinematical ordering (2) for mass sin- 
gularities. However, when dimensional regu- 
larization is used, the implication of the ver- 
tex kinematical ordering gets mixed with the 
fact that the S- and T-channel phase space 
differ by a factor z"^*^. However corrections 
coming from this z-factor do not lead to re- 
ally new structures but are essentially re- 
lated to the anomalous dimensions at lower 
orders. This suggests (see discussion in '*'') 
that these corrections could be an artifact of 
dimensional regularization so that, at the end 
of the calculation, reciprocity is actually re- 
stored leaving RRE unmodified. 

The reciprocity relation (3) can be tested 
for higher order S- and T-anomalous dimen- 
sions 38,39 ^Yiiit one has to accoimt also 
for the arguments of the running coupling 
in S- and T-cases which give contributions 
proportional to the beta-function. However 
beta-function contributions arise also from 
the factorization scheme used for S- and T- 
case and how to account for their reciprocity 
has not been studied yet. Then RRE have 
been tested only for the fixed as contribu- 
tion. In this case RRE can be written as 



7a(Ar) =P(Ar + a7a), 
dz 



(4) 



V{N) 



-f 

Jo 



z^ P{z) 



with 'P{N) a universal function depending on 
as- Eq. (4) has been tested to high order in 
three cases: the non-singlet case, the large x 
and the small x behaviour. 
1) Non-single case. From (4) one has 

^1{N) = V\N) 

+a{jl{N)j\N)+j\N)V\N) 
with 7^ and the n-th expansion coeffi- 
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cients in 0:5 of non-singlet 7cr(-/V) and Va{N). 
Dots are derivatives with respect to N. The 
first is the Gribov-Lipatov relation (indepen- 
dence of a valid only to one loop) . The sec- 
ond is the two loop relation which has been 
pointed out in ^"^ and has been one of the im- 
portant elements used to derive RRE. The 
last one has been verified at three loop or- 
der. 

2) Large x behaviour. The dominant channels 
are the diagonal ones gg and qq. Denoting 
by '^'f\x) the x-space anomalous dimension 
with a = gg or qq, one has 

+4") ln(l-a;)+£»(."^ + --- 

The various coefficients are functions of as- 
The two most singular terms do not depends 
on a. In particular the first term corresponds 
to the classical soft radiation which is uni- 
versal and depends only on the charge A^°''> 
of the source. It can by expressed in terms of 
a physical coupling as A^"'^ = {Ca/TT)a^^^ 
with Ca = Ca,Cf for = gg,qq. Using (4) 
one deduces 

These two relations are verified at three loop 
level, apart for a beta-function contribution 
entering D^"'* which is not considered in (4). 

Very recently the expansion of the lead- 
ing coefficient A has been evaluated in 
SYM with A/'=4 to all order in as (in dimen- 
sional rcgularization scheme suited for super- 
symmetric theories) in the S-case. All coef- 
ficients satisfy the "trascendentality princi- 
ple" (i.e. are given in terms of derivative of 
the Euler psi- function) . The first three terms 
agree with the calculation reported in . Ac- 
cording to (5) the same result should be ob- 
tained for the T-case. 

3) Small X behaviour. The leading order con- 
tributions for N^O (corresponding to small 
x) are given by {as = C^ag/Tr) 



These two expressions satisfy (4). They are 
the result of cancellations in the phase space 
due to coherence of soft radiation: angular 
ordering for the T-casc and transverse mo- 
mentum ordering for the S-case. Thus can- 
cellations in the mass singularity phase space 
are reciprocity related so that RRE incorpo- 
rates them into the universal function P{z). 
RRE can be tested to higher order. One of 
the most interesting outputs is that the "acci- 
dental" absence in the leading BFKL anoma- 
lous dimension of the Og/N"^ and a%/N^ 
terms implies, via reciprocity (4), the fact 
that exact angular ordering is valid beyond 
leading order up to NNLO. 

6 Additional problems in 
hadron-hadron collisions 

Hard processes in hadron collider are initi- 
ated by elementary hard cross sections with 
2 incoming and n outgoing partons 

PXP2 DY, WW, ZZ • • • 

P\P2 1 Pt-Higgs, Higgs-l-jet • • • 

V1P2 — > 2 dijet-distributions, jet-shape ■ • • 

P1P2 n n> 2 , multi-jet distributions 

For n > 2 one has a new QCD challenge. Con- 
sider inclusive distributions with two differ- 
ent hard scales {Q^Qo^ ^qcd) for which 
logarithmic resummations are needed. An 
example is the out-of-cvent-planc energy dis- 
tribution. Because of the complex structure 
of colour matrices for n > 2, soft gluon at 
large angles contribute to single logarithmic 
accuracy. Consider the hard vertex with four 
partons 



Pi P2 P3 Pi 



(6) 



7-(iV) = ^, ^+{N) = h,^m+ms-N) 



Resummation of collinear and infrared loga- 
rithmic contributions coming from radiation 
emitted off the four primary QCD partons 

gives rise to four Sudakov form factors. Each 
factor has the charge of the emitting parton 
and the hard scale identified (as in e+e^) 
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by the single logarithmic analysis. To sin- 
gle log accuracy one has additional contribu- 
tions coming from large angle soft QCD emis- 
sion which arc rcsummcd into a fifth form 
factor^^. To understand why these large an- 
gle soft terms enter only for n > 2 con- 
sider the hard process (6) with matrix colour 
charges T^. Emission of a soft gluon q off 
these four partons is given by the square of 
the eikonal current. Using charge conserva- 
tion {Ti+T2 = T3+T4) one has 

jHl) = {EiTi^,r (7) 
= E, Ti WM + Mq) + Tl Au{q) 

with Tf = Ci the square colour charge and 
Wi{q) the coUinear and infrared divergent 
emission factor leading to the standard Su- 
dakov form factor of the primary parton i. 
The emission factors As and At are infrared 
but not coUinear singular. While are 
proportional to unitV: the exchanged charges 
r2 = (Ti - and = (Ti - T^f are 
colour matrices (6 x 6 in SU{N) for gg — > 
gg). These last two terms, together with 
the corresponding Coulomb phases, give rise 
to the matrix fifth form factor (which needs 
to be diagonalized) . From this discussion it 
is clear that these additional soft, but non 
coUinear, contributions are absent for less 
than four colour particles since there isn't any 
exchanged channel. 

For gg gg there is a puzzle. The rel- 
evant eigenvalues of the soft distribution (7) 
are symmetric under the exchange between 
external and internal space variables: 



ln|--27ri 
In^ 



Another surprising result concerns the con- 
tribution beyond single-log (ncxt-to-next-to- 
leading) to the fifth form factor: it has been 
found that it is proportional to the cor- 
responding one-loop contribution and again 
one reconstructs the physical coupling A'") 
in (5). This result makes possible a variety of 



resummations at next-to-next-to leading or- 
der. 



7 Jets, small-x and all that 

Jets shape distributions. Examples in hadron 
hadron collisions are the out-of-event-plane 
energy distribution and the energy-energy az- 
imiithal correlations. Both reliable predic- 
tions and experimental data ^^'^^ are diffi- 
cult to obtain. The large number of hadrons 
emitted at LHC adds a further algorithmic 
difficulty. A typical kt jet-finder algorithms 
scales as with n the number of parti- 
cles. A good news is that a jet-finder algo- 
rithms was found using techniques devel- 
oped in computational geometry, that scales 
as nlnn. The techniques required to obtain 
reliable jet distribution predictions are well 
imdcrstood: single logarithmic resummation, 
matching with known fix-order NLO results 
and non-perturbative power corrections. Re- 
liable predictions have been obtained when 
no more than four jet are involved, counting 
also incoming jets that is one in DIS and two 
in hardon-hadron collisions. In this last case 
most events have four jets and analytical cal- 
culations become very laborious (see previous 
section) and the use of automate resumma- 
tion becomes essentially imavoidable. An 
additional difficulty in the study of hadron 
emission associated to new physics events is 
that one may want to describe QCD radiation 
in given geometrical regions. In such a case 
one needs to resum non-global logs and 
this requires solving non-linear equations. 
Small-x physics enters a large number of 
processes and gives rise to many interesting 
problems. The connection between BFKL 
and DGLAP equations has been already 
mentioned. Important developments are on 
linear and non-linear small-x regimes. 

In the linear regime there have been ex- 
tensive studies of the higher-order corrections 
which could stabilize the NLO poor perturba- 
tive convergence and construct a framework 
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useful for phcnonicnological applications. 
The general method consists ^'"^ of matching 
small-a; resummation with coUinear singular- 
ity rcsummations (running as, anomalous di- 
mensions, factorization scheme...). 

The non-linear regime, a major issue in 
small- a; physics, is characterized by satura- 
tion which should be seen in heavy ion col- 
lisions A major point is how to go beyond 
the Balitsky-Kovchegov equation viewed as 
a mean field approximation of a more gen- 
eral equation taking into account unitarity 
in perturbative QCD. A number of new for- 
mulations have been proposed, but it seems 
that a complete and consistent formulation 
is still lacking. An important attempt to ac- 
count for general non-linear QCD corrections 
at small-x is based on the observation that 
the BK equation is in the same universality 
class as the Fisher-Kolmogorov-Petrovsky- 
Piscounov equation introduced in statistical 
physics to discuss reaction-diffusion phenom- 
ena. Therefore methods of stochastic physics 
are used ^® to study small-x QCD phenomena 
(and viceversa). 

Perturbative QCD studies of multiple 
interactions in DIS and hadron-hadron colli- 
sions are important for various reasons: they 
are needed to set up the theoretical frame- 
work for the small-x non-linear equations as 
the BK equation and to perform phenomeno- 
logical studies including multi-jet final states 
(background to new physics), heavy flavour 
jets (near forward direction) and underlying 
event 

Diffraction is phenomenologically impor- 
tant since even hard events have difi'ractive 
components. Moreover, diff'raction could be 
used to study Higgs production Diffrac- 
tion is a non-perturbative fields and phe- 
nomenological models are needed. Their 
basis is unitarity and interactions between 
Pomerons 



8 Final remarks 

This (theoretical) review of QCD is incom- 
plete and mostly neglects the phenomenology 
of long distance physics. The list of issues 
here discussed and of the ones not discussed 
(NLO matching with Monte Carlo simula- 
tions, PDF and structure functions, pho- 
ton emission, two photon scattering, prompt 
photons, power corrections, Bjorken and 
Cross-LlewcUyn Smith sum rules, underlying 
events...) speaks for itself of the importance 
of QCD for high energy physics in general. 
QCD for LHC physics poses new problems 
and gives new information on emission such 
as rotation in colour space and consequent 
non-planar corrections. 

Very laborious calculations, needed for 
LHC accurate predictions, reveal simple 
structures and properties. This fact points to 
the possibility of a new more efficient formu- 
lation of QCD. Is it possible that this could 
be provided by ideas developed in string the- 
ory? There are indications in this directions 
although not yet convincing: the abundant 
phenomenological results obtained within the 
framework of AdS/CFT as a way to "solve" 
QCD; Feynman diagrams in twistor formu- 
lations. However, for the moment, all these 
developments are only formal and, to make 
a decisive step in the understanding of QCD, 
one needs to find their physics basis. 
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